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12 Abstract 

13 

14 The accurate quantification of forest 3-D structure is of great importance for studies of the 

15 global carbon cycle and biodiversity. These studies are especially relevant in Africa, where 

16 deforestation rates are high and the lack of background data is great. Mangrove forests are 

17 ecologically significant and it is important to measure mangrove canopy heights and biomass. 

18 The objectives of this study are to estimate: 1. The total area, 2. Canopy height distributions and 

19 3. Aboveground biomass of mangrove forests in Africa. To derive mangrove 3-D structure and 

20 biomass maps, we used a combination of mangrove maps derived from Landsat ETM+, LiDAR 

21 canopy height estimates from ICEsat/GLAS (lee, Cloud, and land Elevation Satellite/Geoscience 

22 Laser Altimeter System) and elevation data from SRTM (Shuttle Radar Topography Mission) for 

23 the African continent. More specifically, we extracted mangrove forest areas on the SRTM DEM 

24 using Landsat based landcover maps. The LiDAR (Light Detection and Ranging) measurements 

25 from the large footprint GLAS sensor were used to derive local estimates of canopy height and 

26 calibrate the Interferometric Synthetic Aperture Radar (InSAR) data from SRTM. We then 

2 7 applied allometric equations relating canopy height to biomass in order to estimate above ground 

28 biomass (AGB) from the canopy height product. The total mangrove area of Africa was estimated 

29 to be 25 960 km 2 with 83 % accuracy. The largest mangrove areas and greatest total biomass was 

30 found in Nigeria covering 8 573 km 2 with 132 xlO 6 Mg AGB. Canopy height across Africa was 

31 estimated with an overall root mean square error of 3.55 m. This error also includes the impact of 

32 using sensors with different resolutions and geolocation error which make comparison between 

33 measurements sensitive to canopy heterogeneities. This study provides the first systematic 

34 estimates of mangrove area, height and biomass in Africa. Our results showed that the 

35 combination of ICEsat/GLAS and SRTM data is well suited for vegetation 3-D mapping on a 

36 continental scale. 

37 

38 1, Introduction 

39 

40 The measurement of forest biomass is crucial for carbon cycle and climate change 

41 studies. However, the amount and distribution of forest biomass is still poorly 

42 understood. Global estimates of terrestrial biomass range from 385 x 10 9 Mg, to 650 x 

43 1 0 9 Mg and forests alone hold about 70-90% of the terrestrial biomass (Houghton et al. 

44 2009). Mangrove forests only cover about 1% of the Earth’s terrestrial surface, but they 

45 are amongst the highest carbon storing and exporting ecosystems globally (Dittmar et al. 

46 2006, Donato et al. 2011). 

47 Estimating distribution and biomass of mangrove forests is challenging due to the 

48 difficult physical environment of these forests. They are constantly inundated by diurnal 
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49 tides and the characteristic aboveground roots often hinder in situ measurements. Large- 

50 scale field measurements of mangroves are therefore rare to inexistent. The 

51 measurements that do exist are usually tailored towards a particular study, and the 

52 sampling and measurement methodologies vary. In Africa, studies of mangroves have 

53 focused on forest composition and zonation (Adams 2004; Dahdouh-Guebas 2004a; De 

54 Boer 2002; Ukpong 1995), management and utilization of mangrove products (Traynor 

55 and Hill 2008; Crona et a/. 2009), the degradation of mangroves (Kruitwagen et al. 

56 2008), and the ecology of mangrove-associated fauna (Faunce and Serafy, 2006). Recent 

57 assessments of mangrove cover in Africa mostly cover small areas, which makes the 

58 comparison with countrywide statistics difficult (Dahdouh-Guebas et al. 2004b). With the 

59 emergence of new remote sensing methodologies, it is now possible to systematically 

60 map mangrove spatial distribution and 3-D structure (Simard et al. 2006, 2008; Fatoyinbo 

61 et al. 2008; Lucas et al. 2007). 

62 Optical Remote Sensing techniques have proven a reliable tool for the estimation of 

63 mangrove forest area globally, as shown by the large number of studies (Aschbacher et 

64 al. 1995; Smith et al. 1998; Dahdouh-Guebas et al. 2000; Kovacs et al. 2001; 

65 Satyanarayana et al. 2001; Dahdouh-Guebas et al. 2002; Sulong et al. 2002; Cohen and 

66 Lara, 2003; Gesche et al. 2004; Wang et al. 2003). The most comprehensive database of 

67 global mangrove cover is maintained by the UNEP World Conservation Monitoring 

68 Center, which published the World Mangrove Atlas (Spalding, 1997). This database is 

69 based on a review of the mangrove literature and mangrove cover estimated from 

70 multiple studies, datasets and methodologies. 

71 For Africa in particular, the data, methodologies and timeframe used to generate 

72 the mangrove maps vary greatly, and a systematic methodology is needed to derive 

73 mangrove cover estimates. An updated version of global maps has recently been 

74 published (Giri et al. 2011). However, to obtain 3D structure and biomass, in addition to 

75 spatial distribution, active remote sensing from LiDAR and InSAR (Interferometric 

76 Synthetic Aperture Radar), are the best measurement tools available. 

77 The only global InSAR and LiDAR datasets currently available are from the 

78 spaceborne SRTM (Shuttle Radar Topography Mission), and ICEsat/GLAS (Ice, Cloud 

79 and land Elevation Satellite/Geosciences Laser Altimetry System). The Shuttle Radar 

80 Topography Mission (SRTM), (Farr et al. 2007), was flown aboard the Space Shuttle 

81 Endeavor in February 2000 (Rodriguez, 2006). The SRTM measured terrain topography 

82 using dual antennae C-band Interferometric Synthetic Aperture Radar (InSAR), covering 

83 areas from 56 0 S’ and 60° N’. SRTM data is freely available at 1-arcsecond (30m) 

84 resolution for the United States and 3-arcsecond (90m), resolution globally. The SRTM 

85 DEM (Digital Elevation Model) is the most accurate, globally consistent elevation dataset 

86 covering 80% of the earth’s landmasses. The SRTM height measurement is in fact biased 

87 by vegetation structure and can therefore be used to estimate canopy height (Kellndorfer 

88 et al. 2004). The GLAS instrument recorded full wavefonn altimetry using a 1064nm 

89 laser that operated from 2003 to 2009. The LiDAR footprint has an approximate diameter 

90 of 70 m, which is separated by 172 m along track (Schutz, 2005). In tropical regions, 

91 sampling is greatly hindered by consistent cloud cover. Although it was primarily a 

92 mission designed for the measurement of icesheet dynamics, it has been used to measure 

93 vegetation structure (Lefsky et al. 2005, 2007; Rosette et al. 2010). Previous work in the 

94 Santa Marta region of Colombia (Simard et al. 2008) has shown the possibility of using 
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95 spaceborne InSAR (Interferometric SAR) and LiDAR data integration, to measure 3-D 

96 vegetation structure and biomass in mangroves. 

97 The objectives of this study are to: 1) estimate mangrove heights on a continental 

98 scale from InSAR and LiDAR integration; 2) estimate the total AGB of mangrove forest 

99 in Africa and 3) estimate the associated errors in our measurements. In this study, we 

100 produce the first continental scale maps of mangrove spatial distribution, 3-D structure 

101 and above ground biomass (AGB), for Africa. We address new challenges introduced by 

102 large-scale mapping that are related to the variety of the biogeographical setting as well 

103 as the accuracy and sampling of data. 

104 

105 2. MATERIALS AND METHODS 

106 

107 2.1. Study areas 

108 In continental Africa, mangroves grow in coastal areas ranging from Mauritania (19° 

109 N’), in the North West to Angola (10° S’), in the South West, and from South Africa (29° 

110 S’), in the South East to Egypt (28° N’), in the North East, including Madagascar. On the 

111 Atlantic Coast of Western Africa there are a total of 7 indigenous species plus one 

112 introduced mangrove palm, Nypa fruticans, which are also found on the Atlantic and 

113 Pacific coasts of the Americas (Spalding et al. 1997). The indigenous species are: 

114 Acrostichum aureum, Avicennia germinans, Conocarpus erectus, Laguncularia 

115 racemosa, Rhizophora harrisonii, Rhizophora mangle and Rhizophora racemosa. The 

116 distribution limit of mangroves coincides with arid regions with rainfall below 30 

117 mm/year (Saenger and Bellan, 1995). 

118 On the Indian Ocean and Red Sea coastlines, the mangrove area is relatively small 

119 compared to the total length of the coastline, due to very arid conditions in areas North of 

120 the equator. There are fourteen species of mangrove present in this area, which differ 

121 from the west coast species. They are: Acrostichum aureum, Avicenna marina , Bruguiera 

122 cylindrical, Bruguiera gymnorrhiza, Ceriops tagai, Excoecaria agaiiocha, Heritiera 

123 littoralis, Lumnitzera racemosa, Pemphis aciduia, Rhizophora mucronata, Rhizophora 

124 racemosa, Sonneratia alba, Sonneratia caseoiaris and Xylocarpus granatum. The largest 

125 diversity on the continent is found in Mozambique, where ten of the species are present 

126 (Spalding et al. 1997). 

127 

128 2.2. Mangrove extent from Landsat 

129 

130 LANDSAT TM GeoCover data was acquired from the University of Maryland 

131 Global Land Cover Facility (http://glcf.umd.edu). The GeoCover dataset consists of 

132 Landsat data that has near global coverage and is available for three time periods ranging 

133 from 1973 to 2001. The Landsat ETM data used in this study had been orthorectified and 

134 georeferenced (Tucker et al. 2004). 

135 A total of 1 17 Landsat ETM+ scenes from 1999 to 2002 were subset to include only 

136 low elevation coastal areas where mangroves may be present. All areas with elevations 

137 lower than 40 m were identified using the SRTM DEM. An unsupervised ISODATA 

138 classification was then applied to each Landsat image subset to discriminate mangroves 

139 from other types of vegetation (Green et al. 1998; Fatoyinbo et al. 2008). The 

140 classification was filtered using previously published maps, the World Mangrove Atlas 
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141 (Spalding et al. 1997), visual inspection and high-resolution imagery from Google Earth 

142 software. The resulting classes were manually combined into a final classification with 4 

143 landcover types (mangrove, other vegetation, bare ground and water). In mangrove 

144 forests in Central Africa, in particular Gabon and the Democratic Republic of Congo, no 

145 cloudless Landsat scenes were available. In these areas with persistent cloud cover, we 

146 had to use cloud free Landsat data from 1989. 

147 There are no local maps with known accuracy or sufficient field data available to 

148 assess relative accuracy. Therefore, we based our estimation of classification accuracy on 

149 an independent and systematic method for selecting validation points. We used points 

150 separated by 900 m (10 pixels), along a North-South running transect. The points were 

151 also spaced by 0.5 degrees longitude for the coast running from Senegal to Nigeria. For 

152 the remaining areas, we used points separated by 900 m along an East-West running 

153 transect, and spaced by 0.5 degrees latitude. We assessed mapping accuracy by visual 

154 interpretation of high-resolution images in Google Earth software. We only used points 

155 that were classified or identified as mangroves on the landcover map or in Google earth. 

156 

157 2.3. Measurement of tree height from LiDAR-InSAR fusion 

158 

159 ICESat/GLAS wavefonns were acquired from the National Snow and Ice Data Center 

160 (NSIDC) website (http://nsidc.org/data/icesat). We used the GLAM (Global land 

161 altimetry) data product to estimate canopy height. A total of 327 waveforms were used to 

162 estimate tree height in this study, as GLAS footprints were not available in all mangrove 

163 areas. GLAS data was available for sites in Senegal, Gambia, Guinea Bissau, Guinea, 

164 Nigeria, Cameroon, Gabon, Congo, Angola, Mozambique, Tanzania, Kenya, Eritrea and 

165 Madagascar. 

166 The GLAM product was produced by fitting up to six Gaussian distributions to the 

167 GLAS LiDAR waveform (Zwally et al. 2003). The shape and position of the Gaussians 

168 distributions describe the canopy vertical structure within the LiDAR footprint. It is 

169 generally assumed that the Gaussian peak furthest from the sensor is the ground return 

170 and the beginning of the waveform signal (i.e. first return with voltage above the noise 

171 level) is the return from the top of the canopy (Harding and Carabajal, 2005). The 

172 cumulative distribution (i.e. percentile) of the energy within the wavefonn is generally 

173 used to describe the vertical distribution of scatterers (e.g. leaves and branches) within 

174 the canopy. The percentile is computed from the beginning of the waveform (i.e. last 

175 return above the noise level). A relative height (rh x ) is defined as the distance between 

176 the point where the percentile energy reaches x and the location of the ground peak 

177 defined as rho (Lefsky et al. 2005, Lefsky et al. 2007). Figure 1 shows an example of a 

178 wavefonn and the location of GLAS footprints used. 

179 We only used data from cloud-free profiles and excluded all waveforms that did not 

180 have suitable data for detennining tree heights. We excluded wavefonns with a single 

181 Gaussian peak, which generally meant the footprint measured water or bare soil areas. 

182 We also excluded wavefonns with low signal to noise ratio (i.e. below 50), which may 

183 have been reflected from clouds, or where Gaussians fits may include noise peaks. We 

184 found high signal to noise ratios up to 300 in the GLAS data. 

185 SRTM version 4 data were downloaded from the Consultative Group for Agricultural 

186 Research (CGIAR). We used 30 SRTM scenes to build a single SRTM DEM covering 
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the coast of Africa mosaic. Using the mangrove landcover map, we masked all non- 
mangrove areas on the SRTM DEM. This resulted in an uncalibrated height map of the 
mangrove areas. In forests, the C-band Radar signal penetrates into the canopy to scatter 
with all forest components and the ground. Thus, the radar height estimate (i.e. radar 
phase center) lies somewhere within the canopy volume, which can be used to estimate 
canopy height (Kellndorfer et al. 2004; Gillespie et al. 2006). Based on the reasonable 
assumption that mangroves are located at sea level, the elevation measured by SRTM (i.e 
phase center) is directly related to canopy height and can be calibrated to estimate the 
canopy height of mangrove forests (Simard et al. 2006). 

The SRTM pixel corresponding to the GLAS shots were extracted (Figure 1). 
Assuming that represents the canopy height, we derived linear regressions between the 
GLAS point’s rh 75 values (relative height of the canopy at the 75 th percentile minus rho), 
and DEM height (Hsrtm) values to determine the regression equation of the form: 

rh 75 = a * HsRTM+b (1) 

Studies of forest biomass worldwide have shown that there is a strong correlation 
between tree size, in terms of diameter and height, and tree biomass. In general, the 
Diameter at Breast Height (DBH), of a tree is the strongest predictor of aboveground 
biomass (Chave et al. 2005). For mangrove forests, a global stand height-biomass 
allometric equation was calculated by Saenger and Snedaker (1993): 

Biomass (Mg ha' 1 ), = 10.8 * Height (m) + 35(2) 

This equation was obtained from 43 field datasets distributed globally (r = 0.59 and 
RMSE= 43.8). To compute total aboveground biomass and aboveground biomass 
distribution of mangroves on the continental scale, we used rh 7 5 and equation (2) to 
derive the biomass values as this equation was computed for a large range of tree heights 
and was derived to be applicable globally. [Insert figure 1 here] 

3. Results and Discussion 

All of the results were calculated and mapped on a per country basis to facilitate 
comparison with previously published results and data distribution. The maps are freely 
available for Google earth software at http://www-radar.ipl.nasa.gov/coastal . 

3.1. Mangrove landcover map 

The total area of mangrove cover in Africa was found to be 25 960 km with 83 % 
accuracy. The five largest mangrove areas were - in decreasing order - in: 1) Nigeria, 2) 
Mozambique, 3) Guinea Bissau, 4) Madagascar and 5) Guinea. The smallest area of 
mangroves is found in Mauritania at 0.4 km 2 . With 8 573 km 2 , Nigeria has the fourth 
largest mangrove area in the world, after Indonesia, Brazil and Australia. The overall 
accuracy of the land cover map was of 83 % considering 10% omissions and 7% 
commissions, based on a total of 540 points (Table 1). The main sources of error in the 
landcover map were due to difficulties distinguishing between mangrove forests and 
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233 other forest types, such as coastal forests or rainforests and the presence of clouds, 

234 especially in the equatorial regions. In Central Africa the map accuracy was much lower, 

235 at 68%, due to the high cloud cover. The landcover maps for Nigeria, Cameroon, 

236 Tanzania and Kenya are presented in figure 2 and the breakdown of mangrove area by 

237 state is presented in table 2. [Insert table 1] 

238 Although it is not our objective to assess changes in the spatial extent of mangroves 

239 over time, it is important to compare our results with previous studies. Overall, the 

240 landcover maps show that the mangrove area of Africa is smaller than the previously 

241 estimated 30 000 km 2 (Spalding et al. 1997, FAO 2007). However the exact estimate of 

242 mangrove area change due to natural and anthropogenic disturbances cannot be 

243 calculated because of the differences in data collection methodologies, the variations in 

244 the definition of mangrove forests and the differences in resolution of the datasets used in 

245 the previous estimates. The large decreases in mangrove area estimates are in part due to 

246 degradation in mangrove area but also due to different definitions of “mangrove areas”. 

247 In many studies, mangrove area was overestimated because it was difficult to 

248 differentiate between mangrove forests and adjacent mudflats, salt marshes, swamp 

249 forests and bare areas using low-resolution data (1 km xl km). The consistent cloud 

250 cover in many tropical areas and poor coverage of optical data. This is the case in many 

251 of the tropical regions, with extreme discrepancies in Congo and Cote d’Ivoire for 

252 example. Furthermore, certain studies include the “mangrove palm” Nypa fruticans as a 

253 mangrove species, whereas other studies do not. In this study we did not include bare 

254 ground and mudflats and also did not count uniform Nypa stands as mangrove areas as 

255 much as was possible. Other very large differences in area measurement such as in Egypt, 

256 Cote d’Ivoire, Sudan, are probably due to a lack of up-to-date studies and remotely 

257 sensed data leading to poor mapping capabilities at the time of the study. 

258 A direct comparison or estimation of the amount and rate of decrease or degradation 

259 in mangrove area throughout Africa is difficult, but we know that mangrove areas have 

260 decreased on the continent due to anthropogenic influences. Over 60% of Nigeria’s 

261 mangrove stands are found in the Niger Delta region, yet studies in the Niger Delta have 

262 shown that mangroves have greatly suffered from the development and rapid increase in 

263 oil and gas exploitation in the area and the resulting pollution by oil spills, rapid 

264 urbanization and dredging of canals, as well as the introduction of the invasive mangrove 

265 palm Nypa fruticans (James, G. K. et al. 2007). In general, decreases in mangrove area in 

266 West Africa are primarily attributed to anthropogenic pressures in coastal regions leading 

267 to conversion of land use for the production of salt and rice, urban and tourism 

268 development, pollution, lack of sustainable resource management and recently, the 

269 development of shrimp aquaculture (FAO, 2007). In eastern Africa, large decreases in 

270 mangrove areas are primarily due to felling for household products and conversion to 

271 urban, agricultural and touristic areas and diversion of freshwater from damming. These 

272 measurement inconsistencies justify the need for a systematic approach to mangrove 

273 mapping as presented in this study. 

274 

275 

276 3.2. Height and biomass measurements 

277 
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278 The GLAS-SRTM calibration regression is shown in Figure 2. The resulting linear fit 

279 between the height estimates from rh 75 and the SRTM DEM is: 

280 

281 rh 75 = 1.07*H S rtm+1.70 

282 

283 The RMSE (Root Mean Square Error) is of 3.55 m. Calibrated canopy height maps 

284 for West and East Africa are presented in Figure 3. In previous studies comparing SRTM 

285 derived canopy height with field and airborne LiDAR data, resulted in RMS errors of 1.6 

286 m and 2.0 m respectively (Fatoyinbo et al. 2008; Simard et al. 2006). Our results are very 

287 similar to these studies. These are the lowest errors that can be achieved using data fusion 

288 of these LiDAR and radar sensors without the incorporation of field validation, [insert 

289 figure 2] 

290 Based on our results, the equatorial areas of Africa are the best suited for the growth 

291 of tall mangroves but not for their expansion, since the actual mangrove area is small in 

292 these countries. Average biomasses per country ranged from 76 Mg ha' 1 in the Republic 

293 of Benin to 178 Mg ha' 1 in Congo. The greatest total biomass values were found in 

294 Nigeria and Guinea Bissau, the lowest in Mauritania. 

295 Previous studies of mangrove canopy height, biomass and distribution have shown 

296 that geographical setting is more important in determining mangrove structure and 

297 distribution than the latitudinal distribution (Fatoyinbo et al. 2008). This is particularly 

298 evident on the African continent, and particularly in West Africa, where a great 

299 proportion of mangroves grow within a small range of latitudes, but the forest area and 

300 structure vary greatly. In Nigeria mangroves are extensive and canopy height can be very 

301 tall, but in adjacent Benin and Togo, their distribution is very limited and canopy height 

302 is short. In Senegal, Gambia, Guinea Bissau and Guinea, mangroves extend very far 

303 inland, up to 160 km in Gambia, but at the same latitudes in East Africa, in Somalia, 

304 Djibouti and Eritrea, mangrove forests are sparse. Estuaries and deltas with extensive 

305 freshwater supply are the most advantageous for mangrove growth, both in terms of 

306 height and extent, and have a much greater influence than latitude. Indeed, all of the 

307 mangrove forests with large areas, tall trees and/or high biomass grow either in estuaries 

308 or in deltas. [Insert figure 3 and table 2] 

309 

310 3.3. Error Analysis 

311 

312 The fact that we used 3 different data sets in this study also increases the incidence of 

313 error in our calculations. In the landcover classification, we observed 83% accuracy, with 

314 17% errors from commissions and omissions from the classification. The systematic error 

315 (i.e. bias) from the calibration equation was low at 1%. 

316 Cloud cover was a major source of error, especially in central African nations, where 

317 cloud cover is persistent. Some systematic but localized errors in the SRTM DEM 

318 resulted in overestimation of tree height and biomass, but also in the omission of 

319 mangrove areas. For example on an island in the Niger Delta, the DEM indicated that 

320 canopy height was 363 m. This is a common error with the SRTM DEM on islands that 

321 may have been caused by difficult SRTM interferometric phase unwrapping (i.e. the 

322 method to retrieve elevation from radar interferometric phase). Because this measurement 
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323 was too high for mangroves this area was omitted from the height and biomass 

324 estimation. 

325 The geolocation error of the GLAS instrument ranges from 4.6 m to 53.4 m 

326 (according to NSIDC), which greatly influences the accuracy of the height measurement, 

327 particularly if the canopy is heterogeneous. The actual height derived from the GLAS 

328 wavefonn may therefore not correspond to the mean canopy height of the SRTM pixel 

329 that is measured. The height estimated from the LiDAR waveform is affected by forest 

330 composition and heterogeneity as canopy shape, reflective properties and the associated 

331 photon interactions all influence the structure of the waveform (Rosette et al. 2010, North 

332 et al. 2010). In addition, the waveform is most sensitive to the footprint center since laser 

333 gain decreases with distance from the center of the footprint. Mangrove forests are 

334 characterized by distinct “zones” that are dependent on the location relative to the coast 

335 or river and that show great heterogeneity in forest structure, type and composition 

336 (Tomlinson, 1994). When the GLAS footprint is close to the border of two zones, this can 

337 result in large discrepancies in height measurement (Figure 4). Although low in species 

338 composition, mangrove forests are very heterogeneous, ranging from tall, dense forests to 

339 very short, sparse and shrubby areas within a few hundred meters. The 70 m GLAS 

340 footprint is not always able to characterize this heterogeneity, resulting in discrepancies 

341 with SRTM measurements. For example when looking at the variance within a seemingly 

342 homogeneous forest in Cameroon, we found that within a single 1 km 2 patch, the canopy 

343 height standard deviation was 5 m, showing that the height within a forest can vary 

344 greatly within a small area (Figure 4). Therefore, since the trees measured by SRTM and 

345 GLAS are not exactly the same, the differences between the height measurements and 

346 what we state as the error of the measurement are inflated. The differences in physical 

347 parameters measured by radar and LiDAR, in addition to differences in resolution also 

348 increase the height and biomass estimation error. These combinations of sources of error 

349 are illustrated in figure 5. [Insert figure 4] 

350 The identification of the ground location within the wavefonn influences the 

351 estimate of the canopy height and therefore also of the biomass. In tidal forests, such as 

352 mangroves, the height of the ground, or of the water level may vary depending on the 

353 tidal level. This may influence the GLAS ground and therefore the relative height 

354 estimates. On the other hand, microtopographic features will most likely average out by 

355 selecting the furthest Gaussian peak as the ground. In this study we chose to use rh 75 as 

356 the height of the canopy as this measurement resulted in the lowest error when comparing 

357 to the SRTM measurement. 

358 The RMS error of equation (3) is of 65.4 Mg ha' 1 . This error is high due to large 

359 variability in the measurements taken and the natural variability of the data set. Since this 

360 is a global equation, it does not take into account local variability in height and biomass. 

361 There is generally a great amount of uncertainty when working with height-biomass 

362 allometric equations. Because height is not the most direct indicator of tree biomass 

363 (Chave et al. 2005), some error is always introduced into the estimate when deriving 

364 biomass from height. To obtain more accurate measurements of biomass from radar and 

365 LiDAR data, it is crucial that more reliable allometric equations be developed as a 

366 function of vertical structure parameters. [Insert figure 5] 

367 

368 4. Conclusions 
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369 

370 Mangroves are one of the most important ecosystems in coastal areas in tenns of 

371 ecology and economy, but they are still being destroyed and degraded at great rates. The 

372 lack of field studies and homogeneous historical data has made the calculation of rates of 

373 change in mangrove cover difficult. In this paper, we produced the first systematic 

374 estimate of mangrove cover, structure and biomass for the entire African continent and 

375 Madagascar. This map can now be used as a baseline as the techniques used in this paper 

376 allows the recalculation and reproduction with updated estimates of canopy height and 

377 allometry in Africa as well as comparison with the rest of the world. 

378 The total area of mangrove forest in Africa for the period of 1999 to 2000 based on 

379 the classification of Landsat ETM+ images is of 25 960 km“ with the largest area found 

380 in Nigeria at 8 573 km 2 and the smallest area in Mauritania with 0.4 km 2 . The overall 

381 accuracy of the map was of 83% considering 10% omissions and 7% commissions. This 

382 overall estimate is lower than previous estimates of mangrove cover in the World 

383 Mangrove Atlas (Spalding, 1997), mostly due to classification errors from high cloud 

384 cover and difficulties in distinguishing between mangroves and adjacent forests. We do 

385 believe that there is an overall decrease in mangrove cover that can be attributed to 

386 deforestation and degradation of mangroves from anthropogenic pressures, however, we 

387 cannot accurately quantify the rate and percent decrease in area because of the 

388 differences in methodology and datasets used between the various published estimates. 

389 Since mangroves are a relatively homogeneous ecosystem that grows on flat terrain at 

390 sea level, the results from this study are some of the most accurate we can expect from a 

391 Radar/LiDAR integration study. The height maps derived from SRTM and GLAS data 

392 confirmed this type of data fusion to measure mangrove canopy height to be appropriate, 

393 with an average RMSE of 3.55 m. This value includes the impact of canopy 

394 heterogeneity on the remote sensing measurement that is not geolocated. Previous studies 

395 using SRTM and LiDAR datasets in Colombia measured canopy height with an accuracy 

396 of 2.7 m (Simard et al. 2008). When similar methods using LiDAR were combined with 

397 field data, the RMSE decreased to 1.6 m in Mozambique (Fatoyinbo et al. 2008). To 

398 achieve even higher accuracy, or lower error, field validation of mangrove height and 

399 biomass calibration should therefore be included in future studies. 

400 Overall, only 327 usable GLAS footprints were found for all mangrove areas in 

401 Africa. This is a very small sample size covering only 0.02 % of the total mangrove area. 

402 This is however, the greatest number of systematic height measurements available. 

403 GLAS was not optimized for vegetation measurement, but as the only spacebome LiDAR 

404 it is the only dataset available for continental-scale studies. We look forward to the future 

405 LiDAR and InSAR missions, which will provide greater coverage over forested areas. 

406 

407 5. Acknowledgements 

408 

409 Dr. Fatoyinbo would like to thank the NASA Postdoctoral Program for funding this 

410 research. The work presented in this paper was conducted at the Jet Propulsion 

411 Laboratory, California Institute of Technology, under contract with the National 

412 Aeronautics and Space Administration NASA, and at the NASA Goddard Space Flight 

413 Center. 

414 


9 



6. References 


415 

416 

417 ADAMS, J.B., COLLOTY, B.M. and BATE, G.C., 2004, The distribution and state of 

418 mangroves along the coast of Transkei, Eastern Cape Province, South Africa. Wetlands 

419 Ecology and Management, 12 , pp. 531-541. 

420 

421 ASCHBACHER, J., OFREN, R., DELSOL, J.-P., SUSELOJ.B., VIBULSRESTH,S. 

422 and CHARUPPAT, T., 1995, An integrated comparative approach to mangrove 

423 vegetation mapping using advanced remote sensing and GIS technologies: preliminary 

424 results. Hydrobiologica, 295 , pp. 285-294. 

425 

426 CHAVE, J., ANDALO, C., BROWN, S., CAIRNS, M.A., CHAMBERS, J.Q., EAMUS, 

427 D., FOLSTER, H., FROMARD, F., HIGUCHI, N., KIRA, T., LESCURE, J.P., 

428 NELSON, B.W., OGAWA, H., PUIG, H., RIERA, B. and YAMAKURA, T., 2005. Tree 

429 allometry and improved estimation of carbon stocks and balance in tropical forests. 

430 Oecologia, 145 , pp. 87-99. 

431 

432 CRONA, B. I., RONNBACK, P., JIDDA WI, N., OCHIEWO, J., MAGHIMBI, S. and 

433 BANDEIRA, S., 2009, Murky water: Analyzing risk perception and stakeholder 

434 vulnerability related to sewage impacts in mangroves of East Africa. Global 

435 Environmental Change, 19 , pp. 227 - 239. 

436 

437 DAHDOUH-GUEBAS, F. DE BONDT, R. ABEYSINGHE, P. D. KAIRO, J. G. 

438 CANNICCI, S. TRIEST and L. KOEDAM, N., 2004a, Comparative study of the disjunct 

439 zonation pattern of the grey mangrove Avicennia marina Forsk. Vierh. in Gazi Bay 

440 Kenya. Bulletin of Marine Science., 74 , pp. 237-252. 

441 

442 DAHDOUH-GUEBAS, F., VAN POTTELBERGH, I., KAIRO, J.G., CANNICCI, S. and 

443 KOEDAM, N., 2004b, Human-impacted mangroves in Gazi Kenya: predicting future 

444 vegetation based on retrospective remote sensing, social surveys, and tree distribution. 

445 Marine Ecology-Progress Series, 272, pp. 77-92. 

446 

447 DAHDOUH-GUEBAS, F., VERHEYDEN, A., DE GENST, W., HETTIARACHCHI, S., 

448 and KOEDAM, N., 2000, Four decade vegetation dynamics in Sri Lankan mangroves as 

449 detected from sequential aerial photography: a case study in Galle. Bulletin of Marine 

450 Science, 672, pp. 741-759. 

451 

452 DAHDOUH-GUEBAS, F., ZETTERSTROM, T., RONNBACK, P., TROELL, M., 

453 WICKRAMASINGHE, A. and KOEDAM, N. 2002, Recent changes in land-use in the 

454 Pambala-Chilaw lagoon complex Sri Lanka, investigated using remote sensing and GIS: 

455 conservation of mangroves vs. development of shrimp farming. Environment, 

456 Development and Sustainability, 42, pp. 185-200. 

457 

458 DE BOER, W.F., 2002, The rise and fall of the mangrove forests in Maputo Bay, 

459 Mozambique. Wetlands Ecology and Management, 10 , pp. 313-322. 

460 


10 



461 DITTMAR, T., HERTKORN, N., KATTNER, G. and LARA, R.J., 2006, Mangroves, a 

462 major source of dissolved organic carbon to the oceans. Global Biogeochemical Cycles, 

463 20, pp.1-7. 

464 

465 DONATO D.C., KAUFFMAN J.B., MURDIYARSO D., KURNIANTO S., STIDHAM 

466 M., and KANNINEN, M., 2011, Mangroves among the most carbon-rich tropical forests 

467 and key in land- use carbon emissions. Nature Geoscience, 4 , pp. 293-297. 

468 

469 FARR, T. G., ROSEN, P.A., CARO, E., CRIPPEN, R. T., DUREN, R., HENSLEY, S., 

470 KOBRICK, M., PALLER, M., RODRIGUEZ, E., ROTH, L., SEAL, D., SHAFFER, S., 

471 SHIMADA, J., UMLAND, J., WERNER, M., OSKIN, M., BURBANK, D., and 

472 ALSDORF, D., 2007, The Shuttle Radar Topography Mission. Review of Geophysics, 45 , 

473 pp. 1-33. 

474 

475 FATOYINBO, T. E., SIMARD, M. WASHINGTON-ALLEN, R. A. and SHUGART, H. 

476 H., 2008, Landscape-scale extent, height, biomass, and carbon estimation of 

477 Mozambique’s mangrove forests with Landsat ETM+ and Shuttle Radar Topography 

478 Mission elevation data. Journal of Geophysical Research, 113 , pp. 1-13. 

479 

480 FAUNCE, C.H. and SERAFY, J.E., 2006, Mangroves as Fish habitat: 50 years of field 

481 studies. Marine Ecology Progress Series, 318, pp. 1-18. 

482 

483 FAO FOOD and AGRICULTURE ORGANIZATION of the UNITED NATIONS, 2007, 

484 Mangroves of Africa 1980-2005: Country Reports, Forest Resources Assessment 

485 Working Paper No 135. Rome, FAO. 

486 

487 GESCHE, K., B. MICHAEL, W. STEFAN, and B. GERALD, 2004, Mapping land cover 

488 and mangrove structures with remote sensing techniques: A contribution to asynoptic 

489 GIS in support of coastal management in north Brazil. Environmental Management, 34 , 

490 pp. 429-440. 

491 

492 GILLESPIE, T.W., ZUTTA, B.R., EARLY, M.K., and SAATCHI, S., 2006, Predicting 

493 and quantifying the structure of tropical dry forests in south Florida and the neotropics 

494 using spaceborne imagery. Global Ecology and Biogeography, 153, pp. 225-236. 

495 

496 GIRI, C. and MUHLHAUSEN, J., 2008, Mangrove Forest Distributions and Dynamics in 

497 Madagascar 1975-2005. Sensors, 8, pp. 2104-2117. 

498 

499 GIRI, C., OCHIENG, E., TIESZEN, L. L., ZHU, Z., SINGH, A., LOVELAND, T., 

500 MASEK, J. and DUKE, N., 2011, Status and distribution of mangrove forests of the 

501 world using earth observation satellite data. Global Ecology and Biogeographv, 20, 

502 pp. 154-159. 

503 

504 GREEN, E.P, CLARK, C.D., MUMBY, P.J., EDWARDS, A.J. and ELLIS, A.C., 1998, 

505 Remote sensing techniques for mangrove mapping. International Journal of Remote 

506 Sensing, 19 , pp. 935-956. 

507 


11 



508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 


HARDING, D. and CARABAJAL, C., 2005, ICESat wavefonn measurements of within- 
footprint topographic relief and vegetation vertical structure. Geophysical Research 
Letters, 32 , pp. 1-4. 


HOUGHTON, R.A., HALL, F.G., and GOETZ, S. J., 2009, Importance of biomass in the 
global carbon cycle. Journal of Geophysical Research, 114 , pp. 1-13 

JAMES, G.K., ADEGOKE, J.O., EKECHUKWU, S., NWILO, P. and AKINYEDE, J., 
2007, Satellite-Based Assessment of the Extent and Changes in the Mangrove Ecosystem 
of the Niger Delta. Marine Geodesy, 30 , pp 249-267. 

KRUITWAGEN, G., PRATAP, H., COVACI, A. AND BONGA, S. W., 2008, Status of 
pollution in mangrove ecosystems along the coast of Tanzania. Marine Pollution 
Bulletin, 56, pp. 1022-1031. 

KELLNDORFER, J., WALKER, W., PIERCE, L., DOBSON, C., FITES, J.A., 
HUNSAKER, C., VONA, J., and CLUTTER, M., 2004, Vegetation height estimation 
from Shuttle Radar Topography Mission and National Elevation Datasets. Remote 
Sensing of Environment, 933, pp. 339-358. 

LEFSKY, M., HARDING, D., KELLER, M., COHEN, W., CARABAJAL, C., DEL 
BOM ESPIRITO-SANTO, F., HUNTER, M., DE OLIVEIRA JR., R. and DE 
CAMARGO, P., 2005, Estimates of forest canopy height and aboveground biomass using 
ICESat. Geophysical Research Letters, 32, pp. 1-4. 

LEFSKY, M.A., KELLER, M., PANG, Y., DE CAMARGO, P. and HUNTER, M.O., 
2007, Revised method for forest canopy height estimation from the Geoscience Laser 
Altimeter System waveforms. Journal of Applied Remote Sensing, 1 , pp. 1-18. 

LUCAS, R. M., MITCHELL, A. L., ROSENQVIST, A., PROISY, C., MELIUS, A. and 
TICEHURST, C., 2007, The potential of L-band SAR for quantifying mangrove 
characteristics and change: case studies from the tropics. Aquatic Conservation: Marine 
and Freshwater Ecosystems, 173 , pp. 245-264. 

NORTH, P.R.J., ROSETTE, J.A.B., SUAREZ, J.C. and LOS, S.O., 2010, A Monte Carlo 
radiative transfer model of satellite wavefonn LiDAR. International Journal of Remote 
Sensing, 31 , pp. 1343-1358. 

RODRIGUEZ, E„ MORRIS, E. and BELZ, J.E., 2006, A global assessment of the SRTM 
performance. Photogrammetric Engineering & Remote Sensing, 723 , pp. 249-260. 

ROSETTE, J.A.B., NORTH, P.R.J., SUAREZ, J.C., LOS, S.O., 2010, Uncertainty within 
satellite LiDAR estimations of vegetation and topography. International Journal of 
Remote Sensing, 31 , pp. 1325-1342. 

SATYANARAYANA, B., MOHAMAD, K.A., IDRIS, I.F., HUSAIN, M.-L., and 


12 


554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 


DAHDOUH-GUEBAS, F., 2011, Assessment of mangrove vegetation based on remote 
sensing and ground-truth measurements at Tumpat, Kelantan Delta, East Coast of 
Peninsular Malaysia. International Journal of Remote Sensing, 32 , pp. 1635-1650 

SAENGER, P. and BELLAN, M.F., 1995. The Mangrove vegetation of the Atlantic coast 
of Africa. (Toulouse: Universite de Toulouse Press). 

SAENGER, P. and SNEDAKER. S.C., 1993, Pantropical trends in mangrove above- 
ground biomass and annual litterfall. Oecologia, 96 , pp. 293-299. 

SCHUTZ, B. E., ZWALLY, H. J., SHUMAN, C. A., HANCOCK, D. and DIMARZIO, J. 
P., 2005, Overview of the ICESat Mission. Geophysical Research Letters, 32, pp. 
L21S01. 

SIMARD, M., ZHANG, K.Q., RIVERA-MONROY, V.H., ROSS, M.S., RUIZ, P.L., 
CASTANEDA-MOYA, E., TWILLEY, R.R. and RODRIGUEZ, E., 2006, Mapping 
Height and Biomass of Mangrove Forests in Everglades National Park with SRTM 
Elevation Data. Photogrammetric Engineering and Remote Sensing, 723, pp. 299-3 1 1 . 

SIMARD, M., RIVERA-MONROY, V.H., MANCERA-PINEDA, J.E., CASTANEDA- 
MOYA, E. AND TWILLEY, R.R., 2008, A systematic method for 3D mapping of 
mangrove forests based on Shuttle Radar Topography Mission elevation data, 
ICESat/GLAS waveforms and field data: application to Cienaga Grande De Santa Marta, 
Colombia. Remote Sensing of the Environment, 112 , pp. 2131-2144. 

SMITH, G. M., T. SPENCER, A. L. MURRAY, and J. R. FRENCH, 1998, Assessing 
seasonal vegetation change in coastal wetlands with airborne remote sensing: an outline 
methodology. Mangroves and Salt Marshes, 2, pp. 15-28. 

SPALDING, M.D., BLASCO, F. and FIELD, C.D. (Eds.), 1997, World Mangrove Atlas . 
(Okinawa, Japan: The International Society for Mangrove Ecosystems). 

TOMLINSON, P.B., 1994, The botany of mangroves . (Cambridge: Cambridge University 
Press) 

TRAYNOR, C. H. and HILL, T., 2008, Mangrove utilisation and implications for 
participatory forest management, South Africa. Conservation and Society, 62 , pp. 109— 
116. 

TUCKER, C.J., GRANT, D. M. and DYKSTRA, J. D., 2004, NASA’s global 
orthorectified landsat data set. Photogrammetric Engineering and Remote Sensing, 703 , 
pp. 313-322. 

UKPONG, I.E., 1995, An ordination study of and leaf size differences in two red 
mangrove swamp communities in West Africa. Vegetatio, 116 , pp. 147-159. 


13 


600 WANG, Y., GREGORY, B., JARUNEE, MICHAEL, T., AMANI, N., JAMES, T., 

601 LYNNE, H., ROBERT, B., and VEDAST, M., 2003, Remote sensing of mangrove 

602 change along the Tanzania coast. Marine Geodesy, 26, pp. 35-48. 

603 

604 ZWALLY, H.J., SCHUTZ, R., BENTLEY, C., BUFTON, J., HERRING, T., MINSTER, 

605 J., SPINHIRNE, J. and THOMAS, R., 2003, GLAS/ICESat L2 Global Land Surface 

606 Altimetry Data V018, 15 October to 18 November 2003. Boulder, CO: National Snow 

607 and Ice Data Center. Digital media. 


14 



